SPIN-LATTICE RELAXATION IN ZERO-MAGNETIC FIELD 
INDUSED BY MOLECULAR REORIENTATIONS 
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The stochastic Liouville method is used to analyze the general problem 
of spin-lattice relaxation in zero-field for molecules undergoing Markovian 
reorientations. 

I. Introduction 

The resolution advantage of zero-field (ZF) studies for orientationatty dis- 
ordered materials is well known. In particular, the novel-pulsed ZF NMR 
and NQR technique [1,2] offers an excellent approach to this problem since 
it removes the orientational anisotropy which produces the broad high-field 
line shapes in solids. In ZF NMR and NOR the signal comes from longitu- 
dinal nuclear magnetization, i.e. the rank-I statisiical tensor. Time-domain 
ZF signals have also been observed using the method of perturbed angu- 
lar correlations of 7-ray cascades [3]. From such experiments it is possible 
to extract information concerning the zero-field spin-lattice relaxation (ZF 
SLR) of rank > 2 statistical tensors. The corresponding response function 
depends on the rates and microscopic details (in the slow-motional regime) 
of molecular reorientations which modulate the antisotropic part of the spin 
Hamiltonian. 

Mathematical techniques have recently been developed which make anal- 
yses of ZF NMR spectra in the complete tumbling regime feasible [4-6]. Our 
purpose here is to extend this theory. We present a general formalism which 
enables us to compute the response function of statistical tensors of arbitrary 
rank-/c irrespective of the models used to describe the Markovian molecular 
reorientations. A compact expression for the correspondtng spectral function 
is obtained, which is valid for the complete tumbling regime. To illustrate the 
use of the theory we calculate the rank— 2 perturbation coefficient of 7-ray 
cascades arising from the quadrupole interaction of a spin— 1 nucleus with an 
axially symmetric electric field gradient. All the calculations in this paper are 
confined to situations that are macroscopically isotropic. There is continuity 
with the formalism of the preceding articles [4-6] and intermediate results 



derived there are assumed to have been looked at by the reader. 



2. Theory 

In ZF the spin Hamiltonian for the problem is 



Here is the ]?-component of a second-rank spin tensor operator in the 
laboratory frame. are components of a ZF splitting tenser expressed 

in the molecular coordinate system(the principal axis frame) and 
are the Wigner rotation matrices describing the transformation between the 
two frames. The explicit form of and will depend on the type of 
interaction. 

For sufficiently large molecules in dense media the stochastic reorienta- 
tional process may be assumed to be Markovian. It then follows that an ap- 
propriate ensemble average spin density operator p{Q,t) obeys the stochastic 
Liouville equation (SLE) [7,8] 

= -iH^pin, t) + Up{^, t) (2) 

at 

where H^p = [H^ p\,fi = 1 and is the stationary Markovian operator 
describing the tumbling process. Eq.(0) must be solved with the initial con- 
dition 

piQ,0) = mm = ^ (3) 

which takes into account the fact that for isotropic systems there is an equi- 
librium distribution of molecular orientations 0(0) = 

The status of the spin ensemble can be discussed in terms of statistical 
tensors state multipole moments ) 

mt)=T.P^'''\^:t)Tkp{I) (4) 

k,p 

where the coefficient p^'^P\Q, t) and the irreducible polarization operator Tkp{I) 
[9] are given by 

/.(^^)(n,t) = rr[/3(n,t)f +(/)], 
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Here Cj^^j^p is a Clebsch-Gordan coefficient. Tfie corresponding response 
and spectral functions, G^^it) and &p{s), are obtained as averages over tfie 
equilibrium (isotropic) distribution: 

G''P{t) = J dVLp^\VL,t), (6) 

^^=^(5) = r G^P{t) exp{-st)dt, (7) 

where the tilde denotes Laplace transformation. 
It thus follows from eqs. (g) and (||),([7|) that 

G'P{t) = Tr{ldnD{n)mDHn)U-^rDlpmn,{I) ^T.Tr[alP{m,{I)] 

q q 

(8) 

where D{Q) is the finite rotation operator. Following refs. [4,10], we multiply 
both sides of eq.(0) by D{n) on the left and by {-l)PD+{n)D^^p{n) on the 
right. In isotropic systems this procedure allows integration over Q in the 
general form as reported in ref. [10]. Through a straightforward extension of 
the derivation described in ref. [10] we obtain a compact differential kinetic 
equation 

al^-kt) = -iH\0)af^\t) - r-' Ul^oiit) - E >*',!<Tf )(t)) , (9) 
where r is the mean lifetime between rotational jumps, 

F^(o) = E(-iri^2;^2-, (10) 

where Q = Q — Q,' (see also refs. [4,10] ). The initial condition for eq. (0) is 
obtained from eqs. (0), (|5) and (@): 

<^f' (0) = (12) 
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Formally eq. (0) is similar to the impact equation which describes gas-phase 
relaxation. Reorientations may be classified as either "weak" or "strong" de- 
pending on the angular jump, with its size set by the function f{Q). The new 
formulation of the problem allows a solution irrespective of this circumstance 
in the general form. 

From that purpose let us re-express o"^^^^ in the form 

^i''\t) = U4'\t)UfKQ{I) (13) 
KQ 

It is easy to see that in this representation the response function can be 
written as 

G^''\t) = u-m<^it'\t)]k,-, (14) 

Then we have, using vector notation, 

X(t) = -(iA+')X(t) (15) 

where the column vector X(t) is constructed from the coefficients o-^'^^\t)]KQ- 

The elements of the evolution, A , and the motivational, " = (1 — P/r , 
operator matrices are [9] 

E {-lf{K,Qr,K-Q/Q'Tr{H{0)fK,QXl))S,,, (16) 

K'Q' 

LMN 

where 

{KiQi; K - QY^'^' = (-l)2^+^'[(-l)Mi+^' - 1][{2K + l){2Ki + l)]^/^x 

k'q' i K Ki K' \ 

= {hiN - 4^Jv), ^5v/!v = / f{^)Dl,^{^)d^ (19) 

To derive eq.([T7|) we have used the Clebsch-Gordan series for the product of 
Wigner matrices [ 9 ]. Eq.([T5|) can be solved by Laplace transformation to 
give 

X(s) = M (s)X(O) (20) 
4 



where M(s) = si + zA It is easy to see from eqs. ([T^ ) and ( [TBI ) that in 
this representation 



X(0) = la'm)]KQ 



'21- 



Eq.(^D|) is particularly suitable for numerical computation of the spectral 
function ([71). The key step in the calculation is the inversion of matrix M: ' 



mi 



k-qi 



2k +1 



Since M has finite dimensions the inversion is readily achieved by standard 
techniques. The result of eq.(21) provides a general recipe for calculating the 
response of the rank— /c statistical tensor on Markovian molecular reorienta- 
tions in ZF. The most severe restriction of the model is that the lattice is 
described only in terms of the orientationat degrees of freedom. 

In the case osotropically rotating molecules, f{Q) = /(cos(/3))/47r^, from 
(|T7|) and (|T9D we obtain 



qqi —Sr^-^n^^' AAA 

KQ KiQi — 2^ ^kq K^q^kqi k~qi^KKiOq-QOqi^Qi 
L 



(23) 



where TqI = VFLnn is the orientational relaxation time of the axial L— rank 



^00 



tensor. To illustrate the use of the theory we consider the case where the 
dominant anisotropic part of the spin Hamiltonian is the axially symmetric 
quadrupote interaction [5]: 



\3 



DQKif2o{I) = 



2 eQV,, 
^ 3 2/(27-1 



-Kif2o{I) 



(24) 



where Q is the nuclear quadrupole moment, 



Ki-- 

As follows from ( [TBI) 



(25) 



^KQ K,Qi - ( 1 



\2I+K 
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I and (23) 

2K,+1)Y/'DqKjC§% 



2 K Ki 

20\ J J J 



[(-i)^+^^-i]V^gQ. 

(26) 
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As can be seen from (22) and (25) the components [(Jf^']K — Q are uncoupled 
from the rest of the vector X and the problem reduces to the inversion of the 
matrix M in " Kq" subspace. It is convenient to calculate &^{s) in the basis 
of eigenfunctions of the operator". In the "i^'g" subspace we have 



-1 



[U Ti]KnKini - Kn Kim " In^KKi^nni: UxqKn- CJ^q K 



(27) 



U is the unitary matrix which makes the submatrix (Yk-qK-q (22) di 



qi 



agonal, 7^ = Tq^^ denotes the eigenvalues. It is easy to see that in this 
representation eqs. (|T^ and (21) give 



1 



Lsi + ik +' 



(28) 



fcO,A;0 



where 



\3 



[2K + l)(2A'i + l)(2n + l)(2iVi + 1] 



c -1 



1/2 



X 









r 








I 






n 





(29) 



Consequently only one element of the inverted matrix (M) is needed to 
calculate the spectral function. 



3. Discussion 

In the fast motional limit, DqT02 <C 1 , taking into account (26)-(28) to 
second order in perturbation theory, we have 

G(^^)(s) = (s + Afc)-V^'"H0), G^'p\t) = exp[-Xkt)p^'p\0) 

where 

3 2 k{k + l)[4I{I+l)-k{k + !)-!] 
h = ^(eQK.) ro2 72(2/ _ 1)2 

in agreement with Abragam and Poundfll] 
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For a simple illustration of the formalism introduced in section 2 we con- 
sider a case with 7 = 1. From (26)- (28) we obtain 

where 

1 16 

Ni{s) = ^Dlis + r,l) + -Dlis + r,^') + (s + r,,'f{s + r,l) 

Eq. ( |30|) is identical to there suh of ref. [5](eq.(24)) which describes the 
ZF NMR spectral function [1,2]. Ordinary NMR coils can only detect rank-1 
tensors. However, by experimental observaiion of 7-ray cascades, it is possible 
to extract information concerning the relaxation of /c > 2 statistical tensors. 
In particular, the measurer anisotropy [3, 11, 12] is proportional to G^^{t). 
Sometimes it is convenient to observe the average correlation of all decays: 

G^0(oo)=r^7f G^O(t)exp(--), 

which is just the Laplace transformation at s = where tn is the mean 
nuclear lifetime. 

The spectral function (30) depends on Te2 and tq^. In the "strong collision" 
model /(il) = and Te2 = = r. In contrast, Tq^ = n{n + 1) A- under 
Debye orientational diffusion with coefficient Dr. Thus in the slow tumbling 
regime the precise form of the angular correlation depends on dynamical 
details of the motion. 

It is interesting to compare our exact result (30) with the approximate an- 
alytical solution of the problem which has been obtained by Linden-Bell [12]. 
It is easy to see (by Laplace transformation of the corresponding expressions 
for G^^(t) which have been given in ref. [12] that in the fast motional regime 
our results coincide. However, her approximation is not sufficient to obtain 
good quantitative agreement wich eq. (30) in the slow-motional limit, when 
DqT02 ~ 1. 
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